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ABSTRACT 

A formalism that simultaneously searches for the monopolar and dipolar 
peculiar velocities is presented. The formalism is applied to (1) the Mark III 
catalogue, (2) Lauer and Postman's Abell cluster catalogue, and (3) Riess et al.'s 
Type la supernova sample. The emphasis is drawn to the monopolar peculiar 
velocities, i.e., peculiar Hubble flows, within these samples. The samples show 
inconsistent pecuhar Hubble flows within a depth of ~ 60h~^ Mpc. Beyond a 
depth of ~ SOh"^ Mpc, the Hubble flows of all samples converge to the global 
Hubble flow to better than 10% at the 2a level. The results are compared with 
theoretical predictions. They at face value disfavor models predicting smaller 
pecuhar velocities such as the tilted Cold Dark Matter model. Limitations of 
the catalogues are discussed, so are ways to improve the catalogues so that 
an accurate map of Hubble flows in our local universe can be drawn and be 
compared with theoretical predictions. 

Subject headings: cosmology: observations — cosmology: theory — Large-scale 
structure of universe 
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1. Introduction 

Peculiar Hubble flows (monopolar deviations from a global Hubble flow) are as 
important as bulk motions (dipolar deviations from a global Hubble flow) in reflecting 
the underlying density fluctuation of the universe. But they have not received the same 
attention as bulk motions have, because the uncertainty of our knowledge about the true 
Hubble constant Hq is often larger than the expected peculiar Hubble flows at large scales. 
One can, however, investigate the variation of Hubble flows within a sample without 
knowing the value of Hq. Moreover, the inferred Hubble flows from samples that extend 
significantly beyond 100h~^ Mpc should be so close to the global Hubble flow that it is 
meaningful to investigate the peculiar Hubble flows within the samples. In doing so, the 
ultimate goal is to map out the variation of the Hubble expansion as a function of depth, 
which is directly related to the underlying density fluctuation (Shi, Widrow and Dursi 1996; 
Turner, Cen and Ostriker 1992), and then to use the peculiar Hubble flows detected to test 
theoretical models. 

Here I first present a formalism that simultaneously calculates the peculiar Hubble 
flows and the bulk motions, along with their errors, within a sample. The formalism is 
applied to (1) the Mark HI catalogue based on the TuUy-Fisher relation as a distance 
indicator (WiUick et al. 1996), (2) The Abell cluster catalogue of Lauer and Postman (1994, 
LP thereafter) with the brightest cluster galaxies (BCGs) as distance indicators, and (3) 
the sample of 20 Type la supernovae (SNe) of Riess et al. (Riess, Press, Kirshner 1996). 
All three samples are deep or fairly deep: the Mark HI catalogue extends to 11,000 km/sec; 
the LP sample is complete up to 15,000 km/sec; the SN sample probes the deepest, 28,000 
km/sec. In other words, Hubble flows defined by these samples should approximate the 
global Hubble flow to a fairly high degree (see table 1 and text). 

The results from the three samples are presented and analyzed in section 3-5. They 
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are compared with theoretical predictions from four representative models (table 1): the 
standard Cold Dark Matter (sCDM) model, a tilted Cold Dark Matter (tCDM) model, a 
vacuum energy A dominated Cold Dark Matter (ACDM) model, and a Cold plus Hot Dark 
Matter (CHDM) model. Finally I conclude in section 6 and comment briefly on improving 
the measurements of peculiar Hubble flows and their potential to test theoretical models. 



2. Formalism 

The peculiar Hubble flow and the bulk motion of a sample can be obtained by 
maximizing the likelihood (Kaiser 1988, 1991) 



1 



{S, - 5Hr, - U,fl 



i\2 



m, 5H) = J{- exp -^^ J ' " , (1) 

where (= r*fi,i = x,y,z) is the position of an object in the sample, Sq its estimated 
line-of-sight peculiar velocity with an uncertainty cTg. U (= UiVi^i = x,y,z) is the bulk 
motion of the sample, and 6H is its peculiar Hubble expansion rate. Maximizing L{Ui, 6H) 
with respect to Ui and 5H gives 
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Sq is related to the real pecuhar velocity field v(r) by 

Sq = Vi{rq)fq + Cq, (5) 

where is a gaussian random variable with an uncertainty cr^. Therefore, 

SH ^SH^"^ + SH^'\ (6) 
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The first term is the noise-free contribution from density fluctuations 



while the second term represents the contribution from noises in data. 
The window function in eq. (|^ is 
W'{r) = 
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Its Fourier transformation is 
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The expected r.m.s. pecuhar Hubble flow for the sample is 
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Similarly, 
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If the true Hubble constant is uncertain and one wishes to investigate the variation of 
the Hubble flow within the sample, it is more appropriate to calculate the r.m.s. deviation 
of the Hubble flow of a subsample (with a rate H2) from the Hubble flow of the whole 
sample (with a rate Hi). Thus 

,(Hi-H2\\ /,SHi + H0-SH2- Ho^2. ,(5Hi^2 ,5H2\2 /(5ifiw(5if2\\ 

"U J +1 Ho ) V Ho Ho Ho ^ ^ ^ Ho ^ V Ho Ho 

(14) 

in which 

((^)(^)) = d%wmhwi-m,\^ (15) 

and 

Subscript "1" refers to quantities of the entire sample, and subscript "2" refers to those of 
the subsample. 



3. Peculiar Hubble Flows in the Mark III catalogue 

The Mark HI catalogue (Willick et al. 1996) compiled the distances and redshifts to 
about 3000 galaxies, based on a template TuUy-Fisher relation from the 36 clusters of Han 
& Mould (HM sample thereafter; Mould et al. 1991; Han and Mould 1992; Mould et al. 
1993). Therefore the underlying Hubble expansion of the Mark HI catalogue is defined by 
the HM sample. Table 1 hsts the expected la deviation of this underlying Hubble expansion 
rate Hi from a global Hubble constant Ho in various models, which is typically 2% to 
3%. (In this and all the following calculations, clusters are assumed to have a gaussian 
random motion with a one-dimension dispersion of 300 km/sec, and all bulk motions are 
in reference to the Cosmic Microwave Background.) Figure 1(a) shows the variation of 
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Hubble flows in the HM sample with 2(T errorbars calculated from eqs. to ([T6|) vs. the 
depth R of subsamples. Each subsample is defined to be the clusters closest to us, and 
R is the distance to the furthest cluster in the subsample. The errors of the variation are 
estimated from the second term of eq. (0). In particular, they include the contribution 
from uncertainties in determining bulk motions. 

The figure shows significant positive variation within 40 to 60/i~^ Mpc, translating 
into significant peculiar Hubble flows H2 — Hq given the smallness of Hi — Hq. In figure 
1(b) the variation is compared with 2(j expectations (with noise included) from theoretical 
models. Clearly the positive (if 2 — Hi) /Hi found exceed the 2a expectation from sCDM 
within ~ 50h~^ Mpc, and exceed the 2a expectations from CHDM, tCDM and ACDM 
within a wider range of subsamples. The variation at scales beyond ~ 60h^^ Mpc, however, 
is perfectly consistent with all models at 2a. Not only the variation of Hubble flows in the 
HM sample is extremely large, its bulk flow is also extremely large, at (740 ± 150 km/sec, 
— 267± 134 km/sec, —360 ± 141 km/sec), compared to the noise-free la sCDM expectation 
of 450 km/sec, or the noise-free la tCDM expectation of 240 km/sec. 

Neither the large variation of Hubble flows within ~ 50h~^ Mpc nor the large bulk 
flows may be surprising. Figure 2(a) projects the positions of the inner 20 clusters (with a 
depth of 56h~^ Mpc) and the direction of their bulk motion on the sky. The four clusters 
at the lower right quadrant of the map, Telescopium at 24/i~^ Mpc, OC3627 at 29h^^ 
Mpc, Pavo II at 37h^^ Mpc and OC3742 at AOh^^ Mpc, all show positive radial peculiar 
velocities. Since they lie roughly between us and the Great Attractor (GA, at ~4200 
km/sec and in the direction of Hydro- Centaurus) (Lynden-Bell et al. 1988), their positive 
radial peculiar velocities may be the result of the gravitational pull of the GA. The clusters 
not in the general direction of the GA are further away from it and thus experience a 
smaller infall into the GA. As a result the subsample shows a large dipole in the direction 
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of the four clusters, plus a large positive monopolar flow. Even the entire HM sample 
with 36 clusters going up to 11,000 km/sec may show the influence of the GA, since from 
figure 2(b), a projection of the outer 16 clusters, no new cluster beyond the distance of 
the GA is added in the same direction of the four clusters, and almost all the outer 16 
clusters are in the opposite hemisphere. To further show the effect of the four clusters, 
figure 3 plots {H2 — Hi) /Hi vs. the depth R of subsamples without the four clusters. The 
positive variation of Hubble flows at 40 to 60h~^ Mpc becomes much less significant and 
no longer exceeds the 2a expectations of the four chosen models. The bulk motion of the 
remaining 32 clusters, however, is still large, at (718 ± 205 km/sec, — 216 ± 140 km/sec, 
—344 ±191 km/sec). Clearly, the HM sample is so sparse and inhomogeneous that spurious 
bulk motions and peculiar Hubble flows can result from structures such as the GA. The 
window functions VF'(k) of the sample, on the other hand, may not fully reflect the density 
waves that give rise to these structures. Therefore the theoretical expectations based on 
these window functions can be too small. This also implies that the assumption that the 
HM sample well represents the true Hubble flow may rather be poor. 

Similar calculations are done to the 277 galaxy groups of Mathewson et al. (1992) plus 
the 11 clusters of Willick (1991) (M+W sample thereafter) in the Mark HI catalogue. The 
Aaronson et al. (1982) sample is too shallow to be included. The Mathewson et al. (1992) 
sample is much denser than the HM sample. It also goes up to ~ 10, 000 km/sec, but 
only covers the southern sky. The Willick (1991) sample is too sparse, but it is not overly 
influenced by the GA (thanks to the large distances of its members from the GA). Figure 
4(a) and 4(b) show the projections of the two samples on the sky. Figure 5(a) shows the 
variation {H2 — Hi) /Hi vs. the depth R of the subsamples of the M+W sample. Figure 
5(b) shows the comparison between this variation and 2a theoretical expectations (with 
noises included). The errors and noise-included theoretical expectations are calculated 
only for two representative points, because both quantities are monotonic, slowly varying 
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functions of R. It should also be pointed out that Hi is the Hubble expansion rate of 
the HM sample because it provides the template TuUy-Fisher relation for the Mark III 
catalogue. Therefore, instead of eq. (pilD, 




The bulk motion of the M+W sample is relatively small, (256 ± 57 km/sec, —343 ± 61 
km/sec, 170 ± 47 km/sec), while the noise-free expectation of sCDM is 670 km/sec, and 
that of tCDM is 361 km/sec. 

There are two conclusions that can be drawn from figure 5: (1) the M+W sample 
shows Hubble flows that are 6% to 8% faster than the Hubble flow of the HM sample, 
with a ~ 4(7 significance; (2) the differences between the HM sample of 36 clusters and the 
combined M+W sample are consistent with theoretical expectations within 2a, although 
the consistency is marginal for models that predict smaller peculiar velocities, such as 
tCDM and ACDM. But because of the poor sampling of the HM sample, whether the first 
conclusion implies a significant positive peculiar Hubble flow in the M+W sample is murky. 

4. Peculiar Hubble Flows in the LP Catalogue 

Lauer and Postman's Abell cluster catalogue is a volume-limited sample that includes 
119 Abell clusters with a redshift of <^ 0.05 and a galactic latitude above 13 degrees on 
both hemispheres. Its advantages over the HM sample are that it is volume-limited, more 
homogeneous, and deeper. As table 1 shows, its Hubble flow is expected to deviate from 
the global Hubble flow by only ~ 1%. The standard candle of the catalogue is taken to be 
the luminosity L of the brightest cluster galaxies as a function of the second parameter a, 
the power index of L as a function of the aperture. The distance to a BCG is taken to be 
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its cosmological redshift [|, i.e., 

r, = {cz, - U,fi - SH^rg)/Ho = {cz, - U,fi)/{Ho + SH,) = {cz^ - U,r\)/H^, (18) 

where Hi is the Hubble expansion rate defined by the entire sample. The estimated 
line-of-sight peculiar velocity is 

where Mq is the absolute magnitude of BCG g, and M^{aq) is the magnitude of the standard 
candle at a = a^. 

When applying eqs. (JlD to @ to the Lauer and Postman's catalogue, because Mq 
and aq depend on and thus Ui and one has to iterate eqs. (P to (H) to obtain a 
self-consistent result. Assuming no peculiar Hubble fiow, Lauer and Postman (1994) got a 
bulk flow of the sample = 477 ± 250 km/sec, Uy = -142 ± 273 km/sec, U, = 635 ± 198 
km/sec relative to the rest frame of Cosmic Microwave Background Radiation. I first repeat 
Lauer and Postman's calculation, log-linearly interpolating table 3 of LP to find out the 
dependence of Mq and aq on r^. I get f/^ = 475 ± 285 km/sec, Uy = —123 ± 308 km/sec, 
Uz = 648 ± 225 km/sec, in very good agreement with LP's result. The errors are estimated 
from eq. (|13|), and verified by Monte Carlo simulations similar to those of LP. 

I then consider the possibility of a non-zero 6Hi = Hi — Hq. Since the distance scale 
is established within the sample itself, 6Hi of the sample is really not a variable but a 
well-defined value. In particular, one can always redefine the true Hubble expansion rate as 
Hq = Hq + 6Hi = Hi so that variation of Hubble fiows in the sample can be investigated 
without referring to Hq, just like in the previous section. 

^The cosmological distortion is small and omitted. See the discussion at the end of the 
section. 
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However, if one simply follows LP's way of calibrating M^,(a), redefining 5Hi to zero 
cannot be done. This is due to the fact that the non-linear relation between Sq and 
M^{aq) — Mq always skews the distribution of Sq to the negative direction when M^,(a) — Mg 
is gaussian as calibrated by regression from the M-a distribution. An unphysical and 
negative 6Hi/Hq will always result for such a calibration, regardless of the value of Hq. For 
the LP sample, 5Hi/Ho ^ -1.2%. 

One way to get rid of the unphysical 6Hi and set 6 Hi to is to calibrate M^..{a) 
by regression from the {Sq/cZq)-a distribution (my calculation shows that the resulting 
Sq/cZq residual is gaussian with a standard deviation of 0.166 at 6% C.L.). The resultant 
bulk motion of the sample is then Ux = 528 ± 285 km/sec, Uy = —272 ± 311 km/sec, 
and Uz = 607 ± 225 km/sec. Uy being the most uncertain component, shows the biggest 
difference from the LP result. The resultant 6Hi converges to zero to a high precision. 

Once M^{a) is calibrated, one can calculate the variation of Hubble flows within the 
sample. Figure 6 shows the variation with 2a error bars vs. the depth of the LP subsamples, 
and its comparison with noise-included 2cr theoretical expectations. The error bars are 
once again calculated using eq. (p^. They are checked with Monte Carlo simulations and 
consistency is found. The figure shows evidence for negative {H2 — Hi)/Hi within a radius 
of ~ 60h~^ Mpc at the 2a level, indicating a negative peculiar Hubble flow at this scale 
given the tiny deviation of Hi from Hq (table 1). But before the negative peculiar Hubble 
flow is explained with overdensities, possible systematic effects have to be considered first, 
as in the case of the Mark HI catalogue. 

Possible biases in the analysis have been discussed extensively in LP (1994). Among 
them are several radially dependent biases that affect the calculation of peculiar Hubble 
flows. The first is the selection bias. A luminosity-limited catalogue introduces an artificial 
Hubble outflow due to the missing low brightness galaxies. This is apparently not the case 
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with LP's catalogue because it is volume-limited. Secondly, the estimated peculiar velocity 
depends on the deceleration parameter go assumed. But for a catalogue extending only 
to z ^ 0.05, the estimated radial peculiar velocity is only changed by ^ 1% if the true qo 
is changed by 0.5. Since the standard candle in LP's catalogue is estabhshed internally, 
influenced mostly by outlying clusters, the effect on the peculiar Hubble flow of inner 
clusters is ^ 1%. A third bias comes from the random peculiar velocities of BCGs due to 
local non-linearities which tend to scatter more BCGs to lower measured redshifts than to 
higher redshifts. But given a typical value of this random radial velocity of 300 km/sec, the 
velocity bias introduced on the 6000 km/sec scale is only 0.5% (LP). Another important 
concern is whether the BCGs of the inner Abell clusters belong statistically to the same 
population of the entire BCG sample. Table 2 hsts the statistical properties of the inner 
BCG subsamples and the entire BCG sample. They are consistent statistically. 

But a problem arises when one tests whether the detection is dominated by a small 
number of clusters. Figure 7 projects the inner 11 clusters (which show the negative 
peculiar Hubble flow at 2a) on the sky. The subsample is very inhomogeneous, due to its 
small size. In particular, there is one (and only one) cluster (A262) that roughly aligns 
with the bulk motion of the subsample. Since the bulk motion and the peculiar Hubble 
flow are not determined independently, the cluster will certainly play a dominate role in 
determining both quantities. If the clusters is excluded, as shown in the small window of 
flgure 6(a), no variation of Hubble flows is found at the 2a level. Therefore, once again, the 
poor sampling of the subsample renders its result susceptible to systematic effects such as 
small scale structures, unrepresented structures, small number statistics, etc. 
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5. Peculieir Hubble Flows in a Type la supernova sample 

Type la supernova samples probe much deeper (~ 10^ Mpc) and have very precise 
distance measurements (~ 5%). Therefore, even for the 20 Type la SNe compiled by Riess 
et al. (1996), table 1 shows that the expected la deviation of its Hubble flow from the 
global one is only ~ 1%. But because of its extremely sparse sampling, its usefulness to 
investigations of peculiar Hubble flows is significantly compromised. Figure 8 shows the 
variation of Hubble flows in the Type la SN sample. No significant variation is found at any 
scale because of its sparse sampling. Since the errors in the plot roughly scale as the inverse 
square root of the number of objects, as more Type la supernovae are being observed with 
good time coverage, they can certainly provide more accurate results at a wide range of 
depth in the future. 

6. Discussion 

After investigating the variation of Hubble flows in the Mark HI catalogue (Willick et 
al. 1996), the Lauer and Postman (1994) sample, and the Type la supernova sample of 
Riess et al. (1996), there are three conclusions that I would hke to draw: 

• Some signiflcant peculiar Hubble flows are found at face value in the Mark HI 
catalogue and at 50 to dOh"^ Mpc in the Lauer and Postman sample. They at 
face value disfavor cosmological models predicting smaller pecuhar velocities, such 
as tCDM. However, because of the sparse and inhomogeneous sampling of the HM 
cluster sample on which the Hubble expansion of the Mark HI catalogue is defined, 
and of the sparse sampling of the Lauer and Postman sample within the QOh"^ Mpc 
scale, the peculiar Hubble flows found are dominated by statistics of a small number 
of clusters, and thus cannot be trusted at their face value. Further supporting this is 
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the fact that the pecuhar Hubble flows found in the different samples are inconsistent: 
the HM sample shows positive pecuhar Hubble flows below ~ 70h~^ Mpc, while the 
M+W sample shows positive peculiar Hubble flows up to ~ 100/i~^ Mpc; the LP 
sample, on the other hand, shows negative pecuhar Hubble flows at ~ 60h~^ Mpc. 

• All samples agree on that the Hubble flows beyond a depth of ~ 80h~^ Mpc are 
true to the global Hubble flow to better than 10% at ~ 2a C.L. In each sample, the 
variation of Hubble flows at a depth of > 80h~^ Mpc, with respect to the Hubble 
flow of entire sample, are found to be < 10% at 2(T. This combined with the small 
expected deviation (1 to 3% at la, see table 1) of the Hubble flow deflned by each 
sample from the global Hubble flow, ensures that the Hubble flows at > SOh~^ Mpc 
conform to the global Hubble flow to better than 10% at ~ 2a C.L. At depth below 
80h~^ Mpc, the limits on peculiar Hubble flows are somewhat weaker. For example, 
at the depth of the Coma cluster {70h~^ Mpc), while the Lauer and Postman sample 
indicates a less than 15% deviation from the true Hubble flow at 2a, all others imply 
a stronger limit-less than 10% deviation from the true Hubble flow at 2a C.L. 

• Since the peculiar Hubble flows detected in the samples are not trustworthy at the 
moment, so is its comparison with theoretical predictions. However, this is not to 
say that any such comparison will be useless, because the major limiting factor to 
a meaningful comparison between models and observations is neither our uncertain 
knowledge of the Hubble constant, nor the uncertainties in the distance scales, but 
rather the imperfect sampling of our local universe. The errors achieved by the current 
Mark HI catalogue are comparable (~ 1% at ~ 60h~^ Mpc) to the differences in 
noisc-frcc model predictions. Therefore, If samples were improved to yield a consistent 
and reliable detections of peculiar Hubble flows, measurements of peculiar Hubble 
flows would certainly be able to test models. 
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While the first and second conclusions represent an attempt to map out Hubble fiows in 
our local universe, which are extremely interesting in their own rights, the third conclusion 
aims at finding a potential tool to gain further information on cosmological parameters, 
and therefore warrants a further discussion. Since errors in the question are roughly 
proportional to the square root of the number of objects in a sample, if a deep (> 100h~^ 
Mpc), homogeneous (covering all directions and depth) and dense (at least 10^ clusters to 
yield a template TuUy-Fisher relation and ~ 10^ clusters and groups within 100/i~^ Mpc to 
be sampled) sample can be assembled with the TuUy-Fisher relation, a reliable detection of 
pecuhar Hubble fiows may be made with errors at the ^ 1% level, sufficient to test many 
models at the ~ 60h~^ Mpc scale. Or, if more than 100 Type la supernovae are adequately 
observed as distance indicators, the errors of their measured peculiar Hubble flows can be 
cut by more than 1/2 from the current level to ~ 1%. Better yet, if the TuUy-Fisher relation 
is calibrated with Type la supernovae, by observations of Type la supernovae in clusters, 
not only the overall Hubble flow of the TuUy-Fisher sample will represent the global Hubble 
flow better, the two samples can be combined to yield better statistics. 
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Table 1. Parameters of four chosen models*and their predicted la deviations from a 
global Hubble flow for the Han & Mould sample (HM), the Lauer & Postman sample (LP), 

and the Type la supernova sample (SN). 



sCDM tCDM ACDM CHDM 



^tot 


1 


1 


1 


1 




1 


1 


0.3 


1 




0.06 


0.06 


0.03 


0.06 













0.2 


h 


0.5 


0.5 


0.7 


0.5 


n 


1 


0.7 


1 


0.9 



(( 



Ho 

Hi — Hp 
Ho 

Hi — Hp 
Ho 



2\ 1/2 



2\ if 2 



SN 



3.1% 


1.8% 


2.1% 


2.4% 


1.2% 


0.7% 


0.9% 


0.9% 


2.0% 


1.0% 


1.1% 


1.3% 



*The power spectra of the models are from Bardeen et al. 
(1986), Bunn and White (1996), Ma (1996) and Sugiyama 
(1995). 
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Table 2. Comparison of statistical properties in the LP 

sample. 



N 


(a) b 


{Mg-M4ag)) ^ 


^ d 

<^Mg-M^{ag) 


Pks ' 


8 


0.61 


0.141 


0.255 


83% 


9 


0.61 


0.113 


0.253 


71% 


10 


0.59 


0.109 


0.239 


84% 


11 


0.59 


0.106 


0.227 


93% 


119 


0.57 





0.245 


39% 



'^The number of clusters in the (sub)saniple. 
''The average a in the (sub)sample. 
'^The average Mg — M*(ag) in the (sub)saniple. 
^The standard deviation of Mq — Mt,{aq). 

^The confidence level that Mq — M^{aq) is gaussian according to a 
K-S test. 
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Fig. 1. — (a) The variation of Hubble flows in the HM sample with 2a error bars, as a 
function of the depth of subsamplcs. (b) The crosses show the variation of Hubble flows in 
the HM sample. Curves are the 2cr expectations (with noises included) from models. Model 
expectations for subsamples are connected for a simpler presentation. 

Fig. 2. — (a) Projection of the inner 20 clusters of the HM sample and their bulk motion 
in galactic coordinates. Solid circles denote clusters with negative radial peculiar velocities, 
open circles denote clusters with positive radial pecuhar velocities. The cross is the direction 
of the bulk motion, (b) Projection of the outer 16 clusters of the HM sample in galactic 
coordinates. The dotted hue is the equator of the sky when the bulk motion of the inner 20 
clusters is chosen as a pole. 

Fig. 3. — The variation of Hubble flows in the HM sample less four clusters in the lower 
right quadrant of figure 2(a), with 2a error bars. 

Fig. 4. — (a) Projection of 277 groups of the Mathewson et al. (1992) sample in galactic 
coordinates, (b) Projection of 11 clusters of the Willick (1991) sample in galactic coordinates. 
Legends are the same as in figure 2. 

Fig. 5. — Similar to figure 1 for the M+W sample. Errors and theoretical expectations are 
calculated only for one subsample and the entire sample. 

Fig. 6. — Similar to figure 1 for the Lauer and Postman (1994) sample. Plot inside the 
window of (a) shows the variation of Hubble flows when cluster A262 is excluded. 

Fig. 7. — Projection of the inner 11 Abell clusters in the Lauer and Postman sample and 
their bulk motion in galactic coordinates. Legends are the same as in flgure 2. 

Fig. 8. — Similar to flgure 1(a) for the sample of 20 Type la supernovae. 
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